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We study thermalization of the early Universe when the inflaton can decay into the Standard Model 
(SM) quarks and gluons, using QCD arguments. We describe the possible formation of the thermal 
plasma of soft gluons and quarks well before the completion of reheating. Relevant interaction rates 
of leading order processes and the corresponding thermalization time scale is presented. We discuss 
hadronization while thermalizing the decay products of the inflaton, when the reheat temperature 
of the Universe is below the QCD phase transition but above the temperature of the Big Bang 
nucleosynthesis. 



I. INTRODUCTION 

After the end of inflation, which is required to solve 
some of the outstanding issues of hot Big Bang cosmology 
[1] , the homogeneous inflaton starts oscillating around its 
minimum. During this period the average oscillations of 
the inflaton mimics the equation of state of a pressureless 
fluid. The inflaton oscillations persist until the inflaton 
decays and fragments in order to provide a bath of rel- 
ativistic species in kinetic and chemical equilibrium. In 
the literature the temperature of the relativistic thermal 
bath acquired from the complete decay products of the 
inflaton is known as the reheat temperature: T r h. The 
equilibration of the bath takes finite time and the whole 
process is known as thermalization of the inflaton decay 
products [2], The time scale of thermalization depends 
sensitively on the inflaton coupling to the decay prod- 
ucts, which is usually not known when the inflaton is 
considered to be a SM and SU{2>) gauge singlet, [3]. 

Note that it is important for Big Bang nucleosynthesis 
(BBN) that at the time when the electroweak interac- 
tion between neutrons and protons freezes out (around 
the temperature of few MeV) the Universe must be dom- 
inated by the SM relativistic degrees of freedom. There- 
fore it becomes important that the inflaton must couple 
to the SM fields if not directly but via non-renormalizable 
interactions *. 

Although a thermal bath of temperature ~ 1 MeV is 
necessary, there is no direct evidence of a thermal his- 
tory beyond the BBN era. Therefore, the reheat tem- 
perature could in principle lie anywhere between rn^ > 
?rh > C(l) MeV, where is the inflaton mass. In su- 
persymmetric (SUSY) theories there is an upper bound 
on reheat temperature arising from thermal production 
of the superpartncr of the graviton: the gravitino [6,7]. 
For the gravitino mass ~ 100 GeV, T rh < 10 9 GeV [7] 



"Thermalization, reheating, and curvature perturbations 
may arise from the known sector of the minimal supersym- 
metric SM (MSSM) flat directions [4], whose coupling to the 
SM quarks and leptons are known, see [5]. 



(assuming that the Universe expands adiabatically from 
T r h onwards). 

The issue of thermalization has been addressed in 
many papers [8-15] and very recently some concrete ideas 
have been put forward in [16,17]. In earlier papers [8,9,13] 
elastic interactions, such as 2 — > 2 scattering and annihi- 
lation have been considered, while in [10,11], it was shown 
that elastic collisions lead to kinetic equilibrium while re- 
distributing the energy density, and 2^3 (particle num- 
ber changing) processes that lead to chemical equilibrium 
take a longer time. In [14], it was suggested that elastic 
scattering followed by prompt decay might lead to rapid 
thermalization. The new idea behind rapid thermaliza- 
tion was proposed in [16], where inelastic scattering such 
as 2 — * 3 process have been invoked from the very be- 
ginning. It was found that thermalization time scale is 
roughly given by the inverse of the inelastic scattering 
rate r^ el . In this paper, we consider thermalization with 
a hadronization scenario, where T r h might be below the 
QCD scale Aqcd ~ 0.2 GeV. Especially when the infla- 
ton scale is as low as ~ 0(1) TeV [18,19], and the 
reheat temperature is below the QCD scale ~ 200 MeV, 
then the importance of thermalization and hadronization 
of the quark gluon plasma (QGP) becomes an interesting 
issue. Further note that the inflaton decay products still 
have energy ranging from 0{m^) down to T r h- The hard 
hitting quarks and gluons must lose their initial momen- 
tum towards the last stages of reheating. Depending on 
whether there is a hadronic bath, or a soft bath of quarks 
and gluons, the inflaton decay products will lose their en- 
ergy differently. Our main goal in this paper is a qualita- 
tive understanding of thermalization and hadronization 
process in a cosmological context with an order of mag- 
nitude estimations. 



II. INFLATON DECAY TO QUARKS AND 
GLUONS 

After inflation has ceased, the inflaton field <f> begins 
to execute coherent oscillations about the minimum of 
its potential. The perturbative decay occurs over several 
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such oscillations t . The initial stage of reheating begins 
right after Hi n f ~ m^, where Hi n f denotes the Hubble 
expansion of the Universe towards the end of inflation, 
and denotes the mass of the inflaton. Note that the 
inflaton oscillations always dominate the energy density 
until the inflaton decay is completed. 




where is the inflaton decay rate, and pi is the energy 
density stored in the inflaton sector, and a is the scale 
factor of expansion. The subscript in denotes the initial 
time which can be taken to be the onset of inflaton oscil- 
lations. The inflaton oscillations dominate until r ~ T^ 1 . 
The inflaton energy goes into a thermal bath of relativis- 
tic particles whose energy density is determined by the 
reheat temperature T r h, given by 

T rh ^. 9 7 1/4 (r M P ) 1/2 , (2) 

where g* is the number of relativistic degrees of freedom. 
The above estimation is simple and it still stands as the 
only valid estimate of the reheat temperature. 

As <p is a gauge singlet with respect to all gauge sym- 
metries under consideration, it follows that for the case of 
QCD with quarks, it can couple only to color singlets, and 
to SU(2) gauge singlets constructed from the left-handed 
quark field which forms a doublet in the standard model. 
In this paper we mainly concentrate upon interactions 

£<t>gg — j^fp^^ v Fp V (3) 

which describes the perturbative decay of the inflaton 
to gluons. Inflatons could also decay into SM quarks 
through 

Clpqq = Jtp~( H ^ qR ' ^ 

where H is the Higgs doublet, is the SU(2) doublet, 
and qn is singlet. The neutral Higgs scalar eventually 
decays into the SM quarks and leptons via Yukawa cou- 
plings. 

We can now estimate the number density of hard hit- 
ting quarks and gluons at time: r > r in , where t,„ is the 
reference time when the inflaton starts oscillating. The 
number density is given by 



^The initial stages of the inflaton oscillations could give rise 
to non-thermal production of bosons and fermions, usually 
known as preheating [20,21]. Alternatively the inflaton con- 
densate could fragment due to running mass [22,23]. It is im- 
portant to highlight that preheating does not give rise to ther- 
malization. Preheating produces multi particle non-thermal 
spectrum at initial stages which still requires to be thermal- 
ized with chemical and kinetic equilibrium. 



n x (r) = 2nt(a in )(l-e- r ^-T , (5) 

where n^{a,i n ) = p c j > (ai n ) /m^, and we collectively des- 
ignate hard hitting quarks and gluons by \. Note that 
we are neglecting here the particle number changing pro- 
cesses in the above estimation. 



III. FORMATION OF THE THERMAL BATH OF 
SOFT QUARKS AND GLUONS 

The inflaton decay rapidly builds up a large number 
density of quarks and gluons. Note that fragmentation 
or hadronization cannot occur in the early stages of the 
inflaton decay, the formation of bound states being pre- 
vented by the IR cutoff imposed by the finite expansion 
rate of the universe. The time scale of expansion during 
inflaton oscillations is t oxp ~ 1/H(t) > which is 

much smaller than the typical time scale of strong inter- 
actions t s ~ GeV" 1 (for » C(l) GeV). However, 
these mechanisms will be important towards the end of 
the inflaton decay, and we will discuss these in section V. 

We now ask the question as to how this dense system of 
hard quarks and gluons will thermalize, what the corre- 
sponding temperature would be and estimate the time to 
reach the final reheat temperature T r h- A rigorous study 
of thermalization would involve solving numerically the 
relativistic Boltzmann equation for the time evolution of 
the quark and gluon densities, with splitting functions in- 
corporated from pQCD, and consequent logarithmic cor- 
rections to the quark and gluon densities. A full calcula- 
tion of the parton cascade evolution needs to be done in 
the expanding space-time background, which is beyond 
the scope of the present work. It is worthwhile at this 
stage to review and recast the arguments for rapid ther- 
malization that appear in references [16,17] in a QCD 
based approach. 

The mechanism therein involves creating a soft ther- 
mal bath of gauge bosons from scattering of hard par- 
ticles accompanied by the emission of a soft gluon from 
one of the legs of the 2^2 diagram. This is the gluon 
Bremsstrahlung. The emitted soft gluon can now be in- 
volved in a scattering process with another quark (or 
gluon), thereby leading to an exponential growth in the 
number density of soft particles. Due to pair annihilation 
processes, the thermal bath in QCD will be composed of 
soft quarks as well as soft gluons. The soft quarks and 
gluons can form a thermal plasma with an instantaneous 
temperature [2] 

S* 1 2 H{r mst )Y lSl Ml\ , (6) 

where H{Ti nst ) > T,p is the Hubble parameter at the 
time when the soft thermal bath is created. This instan- 
taneous temperature reaches its maximum T max < m^, 
soon after the inflaton field starts to oscillate around the 
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minimum of its potential. Also an important point to 
note here is that since the inflaton energy density is still 
dominating p^ w H 2 Mp ~ a -3 / 2 , the instantaneous tem- 
perature falls as 

T inst ~ a - 3 / 8 ( 5 ,- 1/2 ^) 1/4 , (7) 

instead of Ti nst ~ a -1 . Numerical simulations also sup- 
port this argument [24] . In the next section we estimate 
T 

± max- 

For elastic processes, one cannot reduce the average 
energy of the system quickly since no new particles are 
created (we are ignoring the red-shift of the particle mo- 
menta for the moment). The Brcmsstrahlung process 
redistributes the energy among softer constituents and 
these soft particles act as a very efficient sink of energy 
once the bath starts to form. An analogous mechanism 
for heavy-ion collisions was suggested by Shuryak [25] 
that a thermally equilibrated quark-gluon plasma (QGP) 
would be created in a two-step process ("hot-glue sce- 
nario"), wherein the gluonic bath would be created in a 
much shorter time than it took the quarks to thermalize. 

We now proceed to make some quantitative estimates 
for the cross-sections and scattering rates. Note that 
these are order of magnitude estimates, and more accu- 
rate methods such as using Boltzmann equations with 
collision terms taken from perturbative QCD may be de- 
veloped as has been done for the finite size system in 
heavy-ion collisions [26]. However, in the cosmological 
context, it is adequate to first obtain an order of magni- 
tude estimate for a proposed physical scenario. 

The 2^3 processes and higher order (2 — > n) gluon 
branching will lead to the formation of a plasma of soft 
quarks and gluons. When this is eventually formed, ex- 
changed gluons will be screened over a distance scale 
p^ 1 ~ (1/gT) corresponding to the momentum scale 
gT * [27]. In the case of quarks being exchanged in a 
plasma, we can use the cutoff given by the kinetic the- 
ory which gives the quarks an effective mass of O(gT). 
A more consistent perturbative approach in a thermal 
plasma would be to the Braaten-Pisarski re-summation 
scheme [28] , which incorporates effective propagators and 
vertices for the quark and gluons, and amounts to a re- 
ordering of perturbation theory. Since we are interested 
only in order of magnitude estimations here, we choose 
the bare propagators and minimal expressions for the cut- 
off. 

For a sufficiently dilute system of quarks and gluons, 
the concept of a thermal plasma does not make sense, 



t The magneto static gluons are not screened perturbatively, 
and this scale does not apply to them. If magnetic screening 
effect is also taken into account then p^ 1 ~ (l/g 2 T). How- 
ever, for our estimation purposes, a starting point is to take 
the minimal form for the screening of the gluon propagator. 



in which case, we choose to regulate in the infra-red by 
using, at any instant, the inverse of the inter-particle sep- 
aration of the decay products of the inflaton (hard quarks 
and gluons which we denote by \)- This is determined 
by Eq. (5) at any instant of time. For < H(t), we 
expand the right hand side of the expression Eq. (5) in 
terms of T^/H, and with the help of Eq. (6), we immedi- 
ately obtain the number density n x of hard quarks and 
gluons, and from it, the infra-red cutoff 

/ rpi \ -1/3 

It should be noted that the expressions for energy transfer 
will be less divergent than those for the absolute cross- 
section since the energy transfer also tends to zero for 
very soft scattering. This reduces the efficacy of soft 
processes for thermalization. 

As emerged in studies of thermal and chemical equili- 
bration in heavy-ion collisions, inelastic processes (2 — > 
n, where n > 2) are likely to be the dominant mech- 
anism for rapid thermalization and entropy generation. 
Within perturbative QCD itself, it has been shown that 
processes that are higher order in a s can be more im- 
portant for equilibration. Thus, gluon branching from a 
quark or gluon line must be taken into account. 

To achieve thermalization at a temperature well be- 
low ?ri0, one needs to generate many soft particles so 
that while the energy density remains constant, the 
number density rapidly increases. Then the average 
energy per particle also decreases. The creation of a 
soft gluonic bath can proceed via Bremsstrahlung emis- 
sion of a gluon from a quark or gluon participating 
in a scattering process (2^3 process). The en- 
ergy loss due to Bremsstrahlung can be mitigated by 
the Landau-Pomeranchuk-Migdal (LPM) effect [29] in a 
dense medium. This effect has been studied in a QED [30] 
as well as in a QCD plasma, with significant modifica- 
tions in the energy loss profile of a parton jet traversing 
a QCD plasma [31]. If the mean free path of the in- 
cident parton is small compared to the formation time 
of the emitted gluon (which can become large for small 
angle scattering by uncertainty principle arguments) , in- 
terference effects between multiple scattering has to be 
taken into account. This issue has been addressed in cal- 
culations of photon and gluon production rates from an 
equilibrated QGP, as well as in the phenomenon of jet 
quenching. The importance of this effect can be gauged 
as follows. 

Suppose the incident quark of energy E q emits at a 
small angle 9, a gluon {k^ = (u> g ,k z ,k±)) with energy 
u g ~ T inst . The formation time of the emitted gluon 
would be given by the uncertainty principle as [31] 

2w„ 2 

where 9 w k±/ui g . The ratio of the mean free path to the 
formation time is given by 
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(10) 



where cr incl = o\ 



is the inelastic cross-section. 



It has a logarithmic dependence on t from the 
Brcmsstrahlung emission, and a 1/t dependence from the 
exchanged momentum. The extra emitted gluon gives 
one extra power of a s compared to the elastic case lead- 
ing to 



CTinol ~ ~ log — - ) . 

^min \ 



(ii) 



Substituting the above equation into Eq. (10), and using 



the fact that t n 
L 

T 



in St n X 



m<f, 



g 2 T inst log 



T inat/ m ^ we find 

(12) 



If this ratio is much less than one, the LPM effect causes 
suppression of energy loss, and consequently, the hard 
quarks (or gluons) lose their energy at a slower rate. 
This can weaken the functional dependence of dE/dr on 
the incident energy [31], leading to a suppression of en- 
ergy loss, and an increase in thermalization time. From 
Eq. (12), it is evident that for a given 47ra s <C 1, the 
assumption of additive energy loss from successive scat- 
tering applies, since the fraction T inst /m,p is much less 
than unity implying that L/t ^> 1 (the logarithm cannot 
counteract the dependence from the prefactor). 

Continuing with the inelastic process in the limit 
L/t 3> 1 (additive energy loss), the time scale for an 
inflaton decay product to lose an energy ~ by such 
processes turns out to be much smaller than in the elastic 
case, as was shown in [16]. The reason is that the soft 
gluon population (neglecting its annihilation to fermions) 
grows exponentially since each soft gluon produced radia- 
tively also acts as a subsequent scattering center for the 
production of further soft gluons. The bath will be ac- 
tually composed of a certain number of fermions as well, 
due to the forward and backward annihilation processes 
between quarks and gluons. However, for our estima- 
tion purposes, it only matters that hard quarks/gluons 
are scattering off soft particles, whose total number is 
important though not its species content. Due to the 
(almost) exponential increase in the number of soft par- 
ticles, we estimate that the rate of energy loss is given 
by 



dE 
~dr~ 



(13) 



where n g is the number density of the soft quarks (and 
gluons) n g ~ g*T? nst § 



The inelastic scattering rate is then given by 



rincl ~ EdT 



Air 77i0 



(14) 



Note that if we were to consider higher order gluon emis- 
sion (2 or more emitted gluons), Tinoi would be sup- 
pressed by further powers of a s and will not compete 
with the 2-^3 process. Furthermore, emission of gluons 
softer than Ti nstl which would ultimately require han- 
dling collinear divergences arc not considered since those 
ultra soft gluons would have to be re-scattered back up 
to a momentum of Ti nst . As we are interested only in 
the time scale to create a thermalized plasma at roughly 
Ti ns t, this aspect will not concern us here [16]. We will 
return to the point about the 2 — > n processes at the end 
of this section. Continuing with our estimate, the time 
to lose an energy is then 



7"inel 



1 inel 



dE 
dr 



9*a 2 s T? t log 



a 2T L 



(15) 



Note that this result is a factor 47ra s /g* different than 
previous estimates [16], because the lower cutoff is taken 
as the screening mass rather than the typical temperature 
(which would be the energy of the individual constituents 
themselves, and would be ignorant of many-body effects). 
Now we are able to estimate the largest temperature of 
the instantaneous thermal bath before reheating. By us- 
ing H (tinst) - And, Eq. (6), and T = gl /4 T 2 h /M Pl we 
obtain the maximum instantaneous temperature 



trh 



3/4 2 

9* a s 



Mp 



An m<f> 



1/2 



(16) 



Note that T max is couple of magnitudes larger than the 
reheat temperature, but still smaller compared to the 
inflaton mass: T max < m^. 

It is sufficient for us to estimate the thermalization 
time by studying 2^3 processes, since further gluon 
branching (2 — > n, n > 3) will come with additional fac- 
tors of a s <C 1. The lower scattering cross-section for 
such processes renders them higher-order corrections to 
the inelastic scattering rate and the thermalization time 
scale. This conclusion can change if the mean free path of 
the hard particles becomes comparable to the formation 
time of the emitted gluon, i.e., if the LPM suppression is 
severe. 



~ Ti ns t ■ This is not strictly true if one considers the collinear 
behavior of the Bremsstrahlung process, but we can ignore 
log(m^/T) corrections which only speed up the thermaliza- 
tion rate by a numerical factor of order one. 



§ We are assuming, as in [16], that the average energy loss is 
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IV. TOWARDS HADRONIZATION 

Thus far, wc have described a generic scenario towards 
reheating the Universe without bothering too much on 
the reheat temperature itself. An implicit assumption 
was made that the reheat temperature was sufficiently 
above the QCD scale. However, a priori there is no rea- 
son that the Universe could not be reheated at a tem- 
perature lower than the QCD scale and above the BBN 
temperature. Such a scenario is plausible if the inflaton 
scale is sufficiently low. Keeping the reheat temperature 
T r h < Aqcd ~ 200 MeV, we will concentrate upon the 
hadronization issue. 

Naturally if the thermal bath of soft quarks and glu- 
ons is formed near Ti nst ~ T c , where T c is the critical 
temperature (~ 200 MeV), the quarks and gluons will 
rapidly hadronize. However, if the thermal bath is ini- 
tially formed at a much higher temperature than a GeV, 
it will hadronize at a later time when the expansion of 
the universe has red-shifted the energy down to the con- 
finement scale. This is particularly important if thermal- 
ization occurs within a Hubble time H~ x . 

Let us imagine that there already exists a plasma of 
soft quarks and gluons. Since the temperature of the 
plasma is inversely related to the scale factor of the FRW 
metric as T oc a~ 3 / 8 . It follows that during the inflaton 
oscillation dominated phase 



T(t) 
T(r ) 

s(r) 



1/4 



T 3 (t)o 3 (t) 



s(r ) T3(T )a 3 (T ) 



(17) 



where To denotes any reference time. It is interesting to 
mention here that similar relations are derived for the 
evolution of the relativistic plasma of quarks and gluons 
formed in heavy- ion collisions **. 

At time r = r c when T = T c ~ 200 MeV which is 
the critical temperature for hadronization, the fluid is 
dilute enough that strong non-perturbative forces recom- 
binc the quark-gluon soup into colorless hadrons. This 
constitutes a phase transition whose nature depends on 
the number of light quark flavors (from universality ar- 
guments [33]) and the quark masses (which act as ex- 
ternal fields). Assuming the strange quark mass to be 
light on the Aqcd scale, we can say that the hadroniza- 
tion proceeds by a first order phase transition with bub- 
bles of hadronic matter appearing within the quark-gluon 



"There, it is appropriate to adopt Bjorken's idealized hydro- 
dynamical scenario of the quark-gluon plasma as an expand- 
ing relativistic fluid [32]. Approximate longitudinal boost 
invariance and the ideal equation of state for a relativistic 
fluid imply that the temperature and entropy density de- 
pend on the proper time r as T(r)/T(r ) = (to/t) 1 / 3 , and 
s(t)/s(t ) = (t /t). 



plasma, which grows in number till the end of the phase 
transition [34]. The time for this transition to occur is 
denoted by r^; the hadronization time. 

It is possible to estimate Th as follows [35]. Through- 
out the phase transition, the temperature remains con- 
stant at T c , while entropy is decreased in going from the 
quark-gluon to the hadronic phase. This means that 
Eq. (17) cannot be applicable in describing the phase 
transition. It is helpful to imagine the quark gluon 
plasma as a subsystem where the temperature remains 
T(t ) = T{t) = T c ~ 200 MeV for the duration of the 
phase transition. Then, the ratio of the entropy densities 
is simply the ratio of e + p of the plasma, which scales 
as the fourth power of temperature. From Eq. (17), it 
follows that 



s(t) = ejr) +p(r) 

s(r c ) ' /; 



(18) 



In the case of heavy-ion collisions, a rerun of the above 
steps yields the relation s(t)/s(t c ) = (t c /t) 4 ^ 3 . 

This means that the phase transition in the early uni- 
verse occurs at a slower rate than in heavy- ion collisions, 
since the entropy/particle is some fixed quantity in the 
two different phases. We will see that this fact is borne 
out by our final expression for the hadronization time. 
During the transition, the system can be described by 
a mixed phase, with a fraction /(r) in the quark-gluon 
phase (1 — /(t) in the hadronic phase). Since entropy is 
additive, we find using Eq. (18), 



f(r) 



Sqg(T c ) - S h (T c ) 



{f( Tc )s gg (T c ) + [l-f(r c )] 



s h (T c )} x -1 - s h (T c 

T 



(19) 



The entropy densities in the two phases can be traded 
for the respective degeneracies (g qg ~ 37 and gh <~ 3) to 
give 

f(r) = (f(r c )g qg + [1 - f(r c )]g h ^ g h ) . 

g qg - gh v t 



(20) 



At r = Th, the quark-gluon phase is completely 
hadronized. Thus, it follows from the above equation 
that 



Th = 



<kg_ 

9h 



f(r c ) + 1 - /(t c ) 



(21) 



Starting with quark-gluon matter at /(r c ) = 1, we find 
that Th ~ 10r c . In heavy-ion collisions, this relation gives 
Th ~ 6t c tt . 



t1 ln the heavy-ion case, there is a power of 3/4 for the term 
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Coming back to Eq. (18), we observe that r c denotes 
the time for the plasma to cool down to the critical tem- 
perature for hadronization T c . The origin of time is taken 
to be the moment when the innaton begins its oscilla- 
tions. Taking r; n f ~ m^ 1 , we can estimate an upper 
bound on r c as 

r c < T- 1 ~ m- 1 ^ 1 , (22) 

if the universe is to be reheated below the QCD phase 
transition. Here, is the non-renormalizable coupling 
of the innaton to the SM quarks and leptons. For a low 
reheat temperature, a small innaton coupling as set by 
the above equation is therefore necessary. It is known 
that this time in the lab frame must be on the /zs time 
scale. 

Unlike for the case of heavy-ion collisions, where 
hadrons, once formed, decouple chemically and thermally 
from each other after a certain time and stream freely to 
the particle detectors, that may not be the case here. 
Since the hadronization time is short compared to the 
innaton decay lifetime, it follows that the (hard) innaton 
decay products will continue impinging on the newly- 
formed hadrons and cause break-up due to the large en- 
ergy transfer. How efficient this process is depends on 
the number density of these hard particles. One also 
needs some estimate of the rate of energy lost by these 
hard particles traveling now in a hadronic phase rather 
than a soft gluonic bath. It is interesting to compare 
this mode of energy loss for the hard particles to the pre- 
viously considered inelastic scattering process. We also 
expect fragmentation to be important for lowering the av- 
erage energy, since energetic partons can now fragment 
into a large number of less energetic hadrons, which can 
re-scatter inelastically among themselves, or resonance 
decay into lighter hadrons. We comment on the relative 
importance of these processes below. 

V. FINAL STAGES OF HADRONIZATION 

A. Deep-inelastic scattering (DIS) 

The break-up of hadrons by hard-hitting quarks can be 
assessed through DIS. The structure functions at some 
x, Q 2 can be obtained from the parton distributions that 
appear in the DGLAP equation [36] 



where the momentum fraction of the daughter parton 
x = zy, with y the momentum fraction of the parent 
quark. Pij (z) are splitting kernels that could be writ- 
ten down in lowest order pQCD, for example. Rather 
than perform a detailed DGLAP evolution, we wish to 
estimate the time scale for the hard quarks to lose the 
bulk of their energy. From Eq. (5), we can conclude that 
the number density of hard particles at temperature of 
order the hadronization scale (~ 200 MeV) is very low 
(n q < T 4 /™^ ~ 10~ 6 fm~ 3 ), for an innaton mass of 
to^ = 10 5 GcV (relatively small innaton mass which 
gives rise to reheat temperature smaller than Aqcd) and 
the temperature being around the QCD phase transi- 
tion. At the same temperature, the number of hadrons 
is ri/j ~ (m/ l T/ l ) 3 / 2 e~ mh / T . The Boltzmann suppression 
is least for the lightest hadron, namely pions, whose 
mass = 140 MeV is very close to Th, where the 
hadronization temperature Th — 150 MeV. Therefore 
n h ~ Tf^- 1 ~ 10" 1 fm~ 3 . This implies that while 
not many hadrons break up, the hard quarks can lose 
energy quite efficiently due to the large number of par- 
ticipants. The average fractional energy loss of a hard 
quark per hard collision is not much less than 1/2 or 
1/3. The mean free path of the hard quark is approx- 
imately I ~ (jih &dis) ■ Assuming a typical cross- 
section of micro-barns for the integrated deep-inelastic 
cross-section at a momentum transfer Q ~ 1 GcV [37] ^ , 
we find I ~ 10 4 fm. Since these quarks are relativistic, 
this implies that the hard quarks would lose most of their 
energy within a time ~ l/c~ 10~ 19 s. 

B. Fragmentation and hadron- hadron scattering 

The mechanism of energy loss by 2^ 3 inelastic scat- 
tering can still occur, though it is not correct to use one 
gluon exchange as the only interaction when quarks are 
so dilute. If we continue to use the perturbative esti- 
mate, at the temperature Th, from Eq.( 15), we obtain 
Tinci ~ 10~ n s. However, we believe that this would 
not be an accurate estimate because of the strong non- 
perturbative forces active at this separation scale of the 
quarks. In light of these strong non-perturbative forces, 
we may consider the fragmentation process whereby the 
energy of the quark is degraded by conversion to hadrons. 
Late thermalization and hadronization can happen also 
by fragmentation of the hard quarks followed by hadron- 
hadron scattering. For purposes of estimation, we neglect 
the effect of the hadronic medium and using a typical 
proper time of formation (r) ~ 1 — 2 fm/c, the time in 



within square brackets in Eq. (21), which follows from the 
different scaling law of entropy density with time. The frac- 
tional power implies a reduced hadronization time in heavy 
ion collisions. 



^Although, the possibility of larger momentum transfers 
also implies a much shorter time scale, since the exchange 
of high energy gluons (at Q 2 1 GeV 2 ) would include an 
enhancement factor of a 2 /a 2 m [38] . 
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the lab frame (time dilated) is given by [39] 

* fr a g ~<T}^ (24) 

Since the energy of the hard quark E q ~ ~ 10 5 GeV 
and the hadronic mass is about 0.1 GeV, we find if rag ~ 
10~ 18 s. Compared to the time scale for the DIS process, 
this is only slightly slower, and the two processes might 
be comparable within a more detailed treatment. The 
main difference between them is the dependence of the 
time scale on the inflaton mass (logarithmic for DIS, and 
linear for fragmentation). For ~ 10 5 GeV, they are 
almost the same. Since the late hadronization happens 
in a time that is very short compared to the inflaton 
decay lifetime, complete hadronization is achieved only 
when all inflatons have decayed. It should be noted how- 
ever that fragmentation alone is not sufficient since one 
must also take into account that one produces hadrons 
with a large multiplicity and with energies considerably 
above T/j ~ 0.1 GeV. These hadrons will inelastically 
scatter among themselves, and lose energy. Since these 
integrated cross-sections (ahh) are of order millibarns, us- 
ing the formula I = l/(rih<Jhh), it is found that the time 
to reach the equilibrium temperature of the hadron soup 
is of order 10 _22 s. This is much smaller than the time for 
the initial fragmentation, so the overall conclusion is that 
one must compare the fragmentation process against the 
deep inelastic process. For an inflaton mass ~ 10 5 
GeV, this leads to the picture that the late inflaton decay 
products thermalize via an initial parton shower, followed 
by DIS, then fragmentation followed by hadron-hadron 
scattering, within a time much shorter than the inflaton 
decay lifetime. 



(<~ 10/xs) based on the relationship between the maxi- 
mum temperature of the plasma and the reheat temper- 
ature, which should be below Aqcd- 

This state of hadrons may not be immediately stable, 
however, on account of the small cool-down and thermal- 
ization time in comparison to the inflaton decay lifetime. 
Hard inflaton decay products will impinge on this bath of 
hadrons, and cause some break up into quarks and glu- 
ons once again. One can also safely conclude, on account 
of the small hadronization time, that complete formation 
of a hadronic bath will be stable only when all inflatons 
have decayed and the number density of hard particles is 
negligible. The temperature characterizing that state of 
the thermal hadronic bath composed of non-relativistic 
particles can be considered to be the reheat temperature 
T rh . 

In this paper, we have borrowed ideas and expressions 
from relativistic heavy-ion physics to discuss cosmologi- 
cal issues such as the thermalization of the inflaton de- 
cay products. This intermingling is natural especially 
near the hadronization scale where current experiments 
are active, thus the issue of low reheat temperature for a 
QCD plasma could be successively improved by a more 
rigorous treatment of the early non-equilibrium stages 
prior to thermalization. 
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VI. CONCLUSION 



We have studied the perturbative decay of the inflaton 
to quarks and gluons, focusing particularly on the aspect 
of thermalization and energy loss of these hard parti- 
cles in an expanding universe. For a given inflaton mass 
171$, the coupling has to be small enough in order that 
thermalization is achieved before inflaton decay is com- 
pleted. The thermalization time is short, driven prin- 
cipally by 2^ 3 inelastic processes. After the thermal 
plasma is formed, it continues to cool in the background 
of an expanding universe, until a critical temperature 
for hadronization is reached. The laws of expansion arc 
different than those of a relativistic fluid expanding on 
account of it's own pressure, as is the case with heavy-ion 
collisions. 

While we have not examined in detail the hadroniza- 
tion process, which occurs in a mixed phase at roughly 
a constant temperature and is accompanied by a reduc- 
tion in entropy, we have estimated that the hadronization 
time is likely to be on the order of micro-seconds as well 
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